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ABSTRACT 

Context. Transverse loop oscillations and loop contractions are commonly associated with solar flares, but the two types of motion 
have traditionally been regarded as separate phenomena. 

Aims. We present an observation of coronal loops contracting and oscillating following onset of a flare. We aim to explain why both 
behaviours are seen together and why only some of the loops oscillate. 

Methods. A time sequence of SDO/AIA 171 A images is analysed to identify positions of coronal loops following the onset of M6.4 
flare SOL2012-03-09T03:53. We focus on five loops in particular, all of which contract during the flare, with three of them oscillating 
as well. A simple model is then developed for contraction and oscillation of a coronal loop. 

Results. We propose that coronal loop contractions and oscillations can occur in a single response to removal of magnetic energy from 
the corona. Our model reproduces the various types of loop motion observed and explains why the highest loops oscillate during their 
contraction while no oscillation is detected for the shortest contracting loops. The proposed framework suggests that loop motions 
can be used as a diagnostic for the removal of coronal magnetic energy by flares, while rapid decrease of coronal magnetic energy is 
a newly-identified excitation mechanism for transverse loop oscillations. 
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1. Introduction 

This paper explores the relationship between contraction and os¬ 
cillation of coronal loops in connection with solar flares. It is 
motivated by the M6.4 flare SOL2012-03-09T03;53 and the as¬ 
sociated coronal collapse, which have been described in detail 
by |Sim5es et al.| ( [2013| l, hereafter referred to as Paper I. Paper I 
also noted the occurrence of transverse loop oscillations but did 
not explore their relationship to the collapse. 

Hudson| (|2000[) proposed that when a flare or coronal mass 


coronal collapse include contracting magnetic loops (Khan et al. 


2006) |Liu et al.|2009||Liu & Wang|2009l[20T0l [Lmet al.|2012 


Simoes et al. 2013) Shen et al.||2014 i while associated pertu 
bations of the photosphere may be visible as rapid irreversible 
changes to the photospheric magnetic held that lag the start of 
the flare by a few minutes (Sudol & Harvey|2005[ Wang & Liu 


tion events and further stimulates interest in these phenomena 
(e.g. [Aschwanden & Schrijver|2011[ [White & Verwichte|2012[ 


White et al. 2012) Wang et al. 2012) Gosain 2012) Verwichte 
et al.|2013| l. Oscillations of this type are normally interpreted as 
fast kink magnetohydrodynamic (MHD) modes of cylinder-like 
structures (Edwin & Roberts 1983| l and, on this basis, procedures 
have been developed to deduce coronal parameters from wave 
properties (e.g. Nakariakov et al. 1999) Nakariakov & Ofman 
20011. Reviews on MHD waves and coronal seismology are 
given by Nakariakov & Verwichte (2005|l; De Moortel (2005|l; 


Banerjee et al. ( 2007| l; De Moortel & Nakariakov ) 2012 1 . 


ejection (CME) removes energy from a volume of corona, the 
surrounding magnetic held should collapse or “implode” (see 
also jJanse & Low]|2007[|Hudson et al.||2008| ). This can be un¬ 
derstood by noting the equivalence between magnetic energy 
density and magnetic pressure; thus, a reduction in magnetic en¬ 
ergy density creates a partial magnetic vacuum that draws in sur¬ 
rounding plasma and magnetic held. Documented signatures of 


|2010[|Petrie & Sudol|2010[|Johnstone et al.|2012|). 

Separately, transverse oscillations of coronal loops have 
received much attention since their first direct detection 
( [Aschwanden et ak] [1999[ [Nakariakov et al!) [1999[ ). The 
Atmospheric Imaging Assembly (AIA) on board NASA’s Solar 
Dynamics Observatory (SDO) satellite ( [Lemen et al.|p0T^ 
Pesnell et al. |2012|l continues the stream of reported oscilla- 


Oscillations and implosions have been noted together in two 
events other than the one considered here. Using SDO, Sun et al. 
( 2012| l and Gosain ( 2012[ ) noted loop oscillations during and af¬ 
ter a coronal collapse associated with the X2.2 flare SOL2011- 
02-15T01:44, with Gosain ( 2012) 1 estimating oscillation periods 
and applying seismological formulas to estimate coronal Alfven 
speeds and magnetic held strengths. Their main focus, however, 
was the collapse itself and no conjectures were made about the 
relationship between the collapse and the oscillations. A causal 
connection has been suggested by [Liu & Wang ( [2010|l, who ob- 
served active region NOAA 10808 with TRACE ([Handy et al. 
1999] ) at the time of SOL2005-09-08T21:05 (X5.4). They sug¬ 
gested that oscillations could be the result of a contracting loop 
interacting (colliding) with underlying loops, and stated that im¬ 
plosions should therefore be considered as an exciter of trans¬ 
verse loop oscillations. We shall propose a different idea as to 
the origin of the oscillations, which is motivated by the higher 
quality time series available for our event. 

This paper is organised as follows. Section presents the 
observation that motivates this paper. In Sect.j^ we propose that 
loop oscillation and contraction in this event are different aspects 
of a single response to the removal of magnetic energy from the 
corona. Three types of response are identified and this explains 
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why some loops oscillate while others do not. The paper con¬ 
cludes with a discussion in Sect.j^and a summary in Sect.|^ 


2. Motivating event 


On 9 March 2012, the M6.4 class flare SOL2012-03-09T03;53 
occurred in active region NOAA 11429, with flare onset 
recorded at 03:36:45 UT. The flare was noteworthy for the ac¬ 
companying coronal collapse, which was very clearly captured 
by SDO/AIA (Paper I). It was also associated with coronal dim¬ 
ming in the active region, the departure of a coronal mass ejec¬ 
tion and a global EUV wave. 

Two types of oscillations were excited by the flare. Hot 
plasma close to the flare core appeared to undergo compressive 
oscillations with a timescale of approximately 60 seconds. These 
manifested as quasi-periodic pulsations in the EUV, SXR and 
HXR time-series and may have been standing slow waves (see 
Paper I for the observations and more detailed discussion). In ad¬ 
dition, images collected by SDO/AIA at 171 A reveal transverse 
(kink) oscillations with periods between two and five minutes 
in the motions of various coronal loops overlying the flare core. 
These oscillations began when the loops started to contract and 
continued after the contraction. 


The present paper is motivated by the motions of coronal 
loops overlying the flare arcade, all of which contract during the 
flare, and some of which oscillate as well. Figure[2indicates five 
representative loops of interest. To analyse their motion, the ar¬ 
tificial slit shown in Fig. was used to construct a time-distance 
image of intensity, which is presented in Fig. together with 
fitted loop positions. 

The fitted loop positions were obtained by manually iden¬ 
tifying spatial ranges where intensity along the slit is due to a 
single loop plus some slowly varying background, then fitting 
a Gaussian plus a straight line using the Levenburg-Marquardt 
least-squares method ( [Markwardt [2009) 1. The centre of each 
Gaussian was taken as the loop position and the Icr position 
uncertainty was estimated using the covariance matrix, the chi- 
squared measure of goodness of fit and the number of points 
fitted. The typical Icr uncertainty is 20% of the pixel size. 

The shortest, lowest loops (Cl and C2) both contract without 
detectable oscillation, while the three longer, upper loops (LI, 
L2 and L3) oscillate significantly during and after the contrac¬ 
tion. Their periods are 290 s (LI), 190 s (L2) and 150 s (L3). 

The start and duration of the coronal collapse coincide with 
the start and duration of the flare impulsive phase, subject to 
time delays that correspond to outward propagation of a sig¬ 
nal at 300 km/s (see Paper I for the timing analysis). This re¬ 
lationship can be seen in Fig. by comparing the loop motions 
with the bright EUV emission at the base of the slit and with the 
45-100 keV HXR count rate from the Fermi Gamma-ray Burst 
Monitor (GBM, |Meegan et ah 2009 1 (bottom panel of Fig.|^. LI 
has a slow rise phase before collapsing, which may be connected 
to a departing CME, however, this slow rise is significantly less 
apparent for L2 and not discernible in L3. A fourth loop that os¬ 
cillates and contracts is evident just above L3. The motion of this 
fourth loop is very similar to that of nearby L3, which is brighter 
and easier to track. We therefore use L3 as the representative of 
motions in its part of the active region, which is sufficient for 
development of the discussion in this paper. EUV emission in 
the active region makes it difficult to track Cl and C2 reliably 
beyond 1200 s, but that is not a problem since their contraction 
is completed before then. 



X(arcsecs) 


Fig. 1. SDO/AIA 171 A image of active region NOAA 11429 
at 03:03:00 UT on 9 March 2012. Labelled arrows indicated the 
loops studied in this paper and the solid line represents the ar¬ 
tificial slit used to analyse their motions. Intensity is displayed 
using the standard AIA 171 A colour table. 
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Fig. 2. (top) Time-distance image showing contraction and os¬ 
cillation of coronal loops. Intensity has been scaled logarithmi¬ 
cally to highlight the loops and fitted positions are overplotted 
with Icr error bars, (bottom) Fermi GBM count rate in the HXR 
45-100 keV band. 



The questions we wish to answer are what connects contrac¬ 
tion and transverse oscillation for the upper loops, and why do 
some loops oscillate but not others? 
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Fig. 3. Cartoon of the “removal-of-support” mechanism, (a) In 
the initial equilibrium, inward magnetic tension is balanced by 
a supporting outward force from the magnetic pressure gradient, 
(b) A flare below the loop removes magnetic energy from the 
corona, thereby reducing the loop’s support while the tension 
force is initially unchanged, (c) Unbalanced forces accelerate the 
loop and it moves toward a new position where force balance is 
restored. 






Fig. 4. Classes of “collapse and oscillation” that result from re¬ 
moving support from loop in different ways. The black (solid) 
line in each plot shows the loop position as a function of time 
and the red (dashed) line shows the equilibrium position, which 
is moved to drive the system, (a) Impulsively excited loop os¬ 
cillations result if support is removed rapidly compared to the 
oscillation period, (b) Gradual displacement is produced if sup¬ 
port is removed slowly compared to the oscillation period, (c) 
Oscillation during collapse is seen if support is removed on a 
time scale broadly similar to the oscillation period. 


3. Combined model of contraction and oscillation 


We propose that the contraction and oscillation of loops seen in 
Fig.m are two aspects of a single response to removal of mag¬ 
netic energy from the coronal volume below the loops. 

It is widely thought (e.g. |Benz||2008] [Shibata & Magar^ 


2011 1 Su et al. 20131 that the energy required to power flares 


is released from magnetic fields in the corona by rapid mag¬ 
netic reconfiguration, which may be enabled by magnetic re¬ 
connection. The direct coronal magnetic held measurements re¬ 
quired to conclusively prove this are not presently possible, how¬ 
ever, less direct estimates of changes to the coronal free mag¬ 
netic energy using magnetic extrapolations generally support 
this view (|Schrijver et al.|2008[ ling et al.|2010 Sun et al.|2012[ 


[Malanushenko et al.||2014[ |Liu et al.||20i4| l.The hypothesised 

rapid reconfiguration of the coronal magnetic held implies quick 
changes to the j x B Lorentz force acting on the coronal plasma, 
and these changes should have a variety of consequences, in¬ 


cluding coronal collapse and changes to photospheric magnetic 
fields, both of which are observed (references in Sect.[^, as well 
as, we argue below, leading to transverse oscillations of contract¬ 
ing loops. 

Figurej^sketches the sequence producing a coronal collapse. 
Prior to the flare, the loops are in a stable configuration. In this 
initial equilibrium, the magnetic tension (B ■ V)B///o compo¬ 
nent of the Lorentz force pulls magnetic loops inward but this 
is opposed by outward forces that support the loop (Fig.[^). In 
the low-y6 corona, equilibria are well-approximated by magnetic 
force balance, so the dominant force supporting loops against 
contraction is reasonably considered to be the outward magnetic 
pressure gradient -V(B^/2//o) associated with the decrease of 
magnetic field strength with increasing distance from the core 
of the active region. When a flare occurs, energy for the flare is 
removed from the coronal magnetic held and converted to other 
forms e.g. flare radiative emissions. The local decrease in mag¬ 
netic energy corresponds to a decrease in the magnetic pressure 
of the affected volume. At this point, if we assume that a nearby 
loop, not involved in reconnection, has not yet moved and that 
the magnetic field inside it is unchanged, then the magnetic ten¬ 
sion force for this loop is unaffected by the flare. However, the 
magnetic pressure gradient in the loop’s vicinity has changed. 
Consequently, the forces acting on the loop become unbalanced 
(Fig. [^) and the loop must move towards the energy release to 
restore force balance (Fig.[^). These arguments are in important 
respects equivalent to those advanced by Hudson ( 2000) l. Similar 
events will also occur if a CME expels magnetic held from the 
active region, e.g. |Shen et al.| ( |20T^ . 

Since coronal loops have inertia, they do not respond instan¬ 
taneously to changes in their environment. Instead, they behave 
like oscillators. This can be justified physically by noting that 
magnetic restoring forces act to return a displaced coronal loop 
to the shape/position in which it is in equilibrium. These forces 
depend on the instantaneous shape of the loop and on its envi¬ 
ronment but not its velocity, and generally increase in magnitude 
with increasing loop displacements. If the force perturbations 
acted in the opposite sense, i.e. to accelerate the loop away from 
the equilibrium position, then the loop would be unstable, which 
is clearly not the case for the loops identified in Fig. [T] since 
they were observed over several hours and survived perturba¬ 
tion by the flare. When action of restoring forces is combined 
with inertia, standing waves are possible, which make the loop 
displacement at any location oscillate in time. A more detailed 
justification for treating coronal loops as oscillators, based on 
the MHD equations, is given in Appendix [A| 

For damped oscillations about a fixed equilibrium position, 
a simple harmonic oscillator obeys 


(fix 


dx 


—TT + Cl) {X - Xq) + 2cl>K - = 0, 

dfi df 


( 1 ) 


where x represents displacement of the oscillator, xq is a constant 
specifying the equilibrium position, t is time, w is the frequency 
of the corresponding undamped oscillator and k is the damping 
ratio (a- < 1 for underdamping, /c = 1 for critical damping and 
A- > 1 for overdamping). The same equation can also be rigor¬ 
ously derived for standing kink oscillations of particular equi¬ 
libria as outlined in Appendix ^ Here, the damping term is in¬ 
cluded to represent decay mechanisms such as resonant absorp¬ 
tion ( |Ruderman & Roberts|2002l|Goossens et al.|2002] l: it does 
not affect the conclusions of this paper, but we include it because 
the observed oscillations (Fig. clearly do decay. We also note 
that the frequency of oscillation for curved coronal loops de¬ 
pends (weakly) on the polarisation of the oscillation with respect 
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to the loop’s curvature (|Van Doorsselaere et al.||2004| |Wang &| 
Solanki|2004| , so CO should be set according to the polarisation. 

The purpose of this paper is to show that time-dependent 
changes to the equilibrium position of an oscillator lead to the 
precisely the observed oscillation and displacement behaviours 
seen in Fig.|^ We therefore modify Eq. Q by making xq a func¬ 
tion of time, xo(f)- This modification is introduced in an ad hoc 
manner rather than by derivation from the full MHD equations, 
but it will be justified a posteriori by comparison of the resulting 
solutions to the observed loop motions. This means that we will 
use the equation 


d^x 


dx 


^ -H CO (x - Xo(f)) + - 0, 


df- 


dt 


( 2 ) 


as a guide to the expected coronal loop motions, emphasising 
that real loops and this model equation share the key physical 
properties of inertia, oscillation and changing equilibrium. 

Three different types of response are found depending on 
how the loop’s period of oscillation compares to the time scale 
over which support is removed. These are shown in Fig. 
Driving functions and parameters that produce examples of each 
behaviour can be found in Appendix [B| 

Figure illustrates the type of motion that results if sup¬ 
port is removed rapidly compared to the period of oscillation. 
The solid black line shows the loop position and the dashed red 
line shows the equilibrium position. At some time, the equilib¬ 
rium steps to a new value and the loop accelerates away from 
its original location. On reaching the new equilibrium position, 
the loop overshoots, and thereafter oscillates about the new equi¬ 
librium. After the initial excitation, the amplitude decays due to 
physical damping or other decay processes such as resonant ab¬ 
sorption. In the model equation, making xq a step function leads 
to a straightforward analytic solution where x is a cosine multi¬ 
plied by an exponentially decaying envelope. This scenario can 
be labelled as an impulsively excited loop oscillation. 

The opposite extreme is illustrated in Fig.|^. Here, the equi¬ 
librium changes on a time scale that is much longer than the 
period of oscillation. In this case, the loop’s position is always 
close to the equilibrium position, although lagging it slightly. 
The outcome is that the loop effectively passes through a series 
of equilibria and no oscillation would be detected. This can be 
referred to as gradual displacement. 

The final type of response applies to the intermediate case in 
which the equilibrium position changes on a time scale broadly 
similar to the period of oscillation (i.e. not differing by more than 
an order of magnitude). An example of this is shown in Fig. B- 
When the equilibrium position starts to move, the loop, having 
a certain inertia, requires time to accelerate before it can follow 
the equilibrium position, which changes as rapidly as dictated 
by removal of magnetic energy/pressure. The loop continues to 
accelerate until it overtakes the equilibrium position and the di¬ 
rection of the acceleration is reversed. Thus, the loop oscillates 
as in the impulsive case (Fig. |^) but this time the oscillation is 
superimposed with the collapse/displacement that excited it. 

Comparing the observations shown in Fig.|^to the responses 
illustrated in Fig. the lower lying loops (Cl and C2) fit the 
pattern of “gradual displacement” sketched in Fig. |^. The mo¬ 
tion of L3 could serve as an archetype of “oscillation during col¬ 
lapse” scenario of Fig. and L2 also belongs to this category. 
Finally, LI appears to be at the borderline between “oscillation 
during collapse” and the “impulsively excited loop oscillation” 
scenario of Fig. |^. A more distinct example of “impulsively 
excited loop oscillation” associated with coronal collapse may 
be found in the uppermost loop studied by |Gosain| (|2012|l for 


SOL2011-02-15. The ordering of the different types of motion 
also agrees with the model since loop period increases with loop 
length from C2 to L1. 

Experiments with the model equation have shown the great¬ 
est amplitude oscillations are obtained when the change in equi¬ 
librium position is initially sharp, thereby allowing the equilib¬ 
rium position to rapidly pull away from the loop and subjecting 
the loop to a greater acceleration than if the equilibrium position 
changed more smoothly. Eurther detail is given in Appendix ^ 
This indicates that oscillations are more likely to be detected if 
the process removing magnetic energy from the corona (presum¬ 
ably magnetic reconnection) switches on over a time scale much 
shorter than the period of oscillation, as oppose to ramping up 
over several oscillation periods. 


4. Discussion 


Loop motions of the type described here are interesting as some 
of the most direct evidence for rapid reduction of coronal mag¬ 
netic energy during flares. They also offer some interesting diag¬ 
nostics. 

Eirst of all, our model gives initial and long term displace¬ 
ments toward the coronal volume from which magnetic energy 
has been removed. The loops studied in this paper move toward 
the underlying flare arcade, suggesting that the energy fuelling 
the flare had been stored in the low corona close to the arcade 
( |Sun et al.|2012| also see). This is reasonable since the core of 
active region is where the coronal magnetic fields are strongest 
and therefore can most easily build up significant free magnetic 
energy. It is also consistent with the observation that changes to 
the non-radial photospheric magnetic held at the time of flares, 
which indicate a change in coronal free magnetic energy, are 
usually concentrated near the polarity inversion line, also sug¬ 
gesting that flares are typically powered by free magnetic energy 
stored in the low corona. 

Our model allows for loop motions other than vertical; in 
particular, if a flare were to remove energy from a volume at the 
same height as the non-flaring magnetic loops, then the loops 
would move perpendicular to their plane as force imbalance ac¬ 
celerates them towards the flare volume. That would excite os¬ 
cillations perpendicular to each loop’s plane rather than within 
it (see Wang & Solanki ( 2004| l for discussion of differences be¬ 
tween in-plane and out-of-plane kink oscillations). [White et al.| 
( 2013| l recently presented an interesting event that fits this pat¬ 
tern. Following an Ml.4 class flare, two coronal loops were seen 
to oscillate, initially moving toward one another. A net inward 
displacement of the loops also occurred during their oscillation. 
This behaviour is consistent with our model if magnetic en¬ 
ergy was removed from the volume between these loops. Indeed, 
post-flare loops for this event do form between the eastern foot- 
points of the much longer oscillating loops, in keeping with such 
a picture. 

As well as indicating volumes where magnetic energy has 
decreased, loop motions provide information on how rapidly 
magnetic energy was removed and how rapidly the energy re¬ 
lease switched on. Similar information can already be gained 
from HXR, EUV and radio emissions, however loop motions are 
potentially valuable as an additional diagnostic since they reflect 
changes to coronal magnetic energy without assumptions about 
conversion to other forms. For the event considered in this pa¬ 
per, inspection of Fig.|^shows that most of the collapse occurred 
within about 300 s, with the entire contraction completed within 
900 s. The logical inference is that the removal of magnetic en¬ 
ergy also followed this pattern. These values are consistent with 
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what one would deduce from the Fermi GBM 45-100 keV HXR 
emission (bottom panel of Fig. |^. Meanwhile, the large ampli¬ 
tude of the oscillations suggests that the switch-on time for the 
energy release is at most a quarter of the shortest observed pe¬ 
riod, i.e. 40 sec (see Appendix[C]for analysis). This time scale is 
comparable to the shorter rise times of the HXR emission. 

The findings of this paper are also of interest to the study 
of coronal oscillations and seismology. In particular, removal of 
magnetic energy from near a coronal loop joins other mecha¬ 
nisms for excitation of kink oscillations, such as interaction with 
a fast wave shock (blast) wave ( |Nakariakov et al.|1999[ |Pascoe| 
jet al.|200^ , reconnection outflow ( White et al.|2012|l, or vortex- 
shedding in response to a flow relative to the loop ([Nakariakov 


et al.||20()9]l. For the event in this paper, the other mechanisms 


can be ruled out because the loops initially move toward the 
flare rather than away from it (ruling out the blast mechanism), 
the loops exist before the flare/eruption and remain cool (in con¬ 
trast to post-reconnection loops), and the oscillations have their 
maximum amplitude at the start of the collapse (indicating that 
vortex-shedding is not responsible here). Comparative motion of 
different loops also indicates a driving signal moving away from 
the flare at around 300 km/s, which is consistent with informa¬ 
tion about the removal of magnetic energy propagating out from 
the flare site at the fast speed (for low-/l plasma the fast speed 
is approximately the Alfven speed) but seems inconsistent with 
driving from above, e.g. by the departing CME. 

In what situations will removal-of-support excite oscillations 
and how common do we therefore expect these events to be? 
Firstly, the mechanism requires that magnetic energy is con¬ 
verted on a time scale broadly comparable to or shorter than 
the oscillation period of the loop. Taking the impulsive phase 
of a flare as a proxy for removal of coronal magnetic energy, 
flare impulsive phases typically last between several tens of sec¬ 
onds and tens of minutes ( [Fletcher et al.|20lT| l. Meanwhile, ex¬ 
amples of transverse loop oscillations collected by jAschwandenj 
et al. I (2002 1 and White & Verwichte (20121 show periods in the 


range 1.7 to 33 minutes (the lower limit is probably observa¬ 
tional since SDO/AIA has a mean cadence of 12 seconds, ex¬ 
cellent by historical standards, and a minimum of six points per 
period is required to convincingly resolve an oscillation). In gen¬ 
eral terms, then, the typical duration of the flare impulsive phase 
is comparable to the lower end of loop periods reported in the 
observational literature, which suggests coronal collapse should 
commonly excite oscillations. A significant proportion of these 
will be of the “oscillate during collapse” type while the longest 
period loops will better fit the “impulsively excited” template. A 
second consideration is the amount of magnetic energy removed 
from the corona. This paper has referred to events with flares of 
class Ml.4 ([White et al.|20T3] l, M6.4 (Sect.[2|, X2. 2 ( [Sun et al.[ 
[2012[ [Gosain[[2012[ l and X5.4 ( [Liu & Wang 2010[ ), so a broad 
range of M and X flares seem to excite oscillations in this way, 
the main difference being that larger flares perturb more of their 
active region and produce larger loop displacements. Finally, we 
have demonstrated that large amplitude oscillations are favoured 
by rapid switch-on of the energy release. This may be the limit¬ 
ing factor, hence we recommend that anyone seeking new exam¬ 
ples first look to flares where the impulsive phase commenced 
suddenly. 

Since many loops in the same active region may be excited 
by coronal collapse (as in the SOL2012-03-09 and SOL2011- 
02-15 events), oscillations excited by coronal collapse could be 
particularly useful for multi-loop coronal seismology, e.g. us¬ 
ing oscillations of several loops to build a picture of the coro¬ 
nal Alfven speed throughout an active region. This application 


is complicated slightly by the need to separate the loop motion 
into oscillation and displacement components, however, we have 
found this issue surmountable, as will be shown in future work. 

Finally, the conclusions of this paper rely only on coronal 
loops having inertia, being natural oscillators and having equi¬ 
librium shapes/positions that change during flares. These points 
are well established observationally and no further assumptions 
are required, so we expect the simple analysis in this paper to be 
a good guide to the qualitative behaviour. Nonetheless, it could 
still be informative to perform 3D MHD simulations of the loss- 
of-support mechanism, which would provide further validation 
and add details about the precise manner in which the loop’s 
equilibrium shape/position changes in response to the removal 
of coronal magnetic energy. We anticipate such simulations in 
the near future. 


5. Conclusions 

This paper has argued that coronal loop contraction and oscilla¬ 
tion in flares can occur as part of a single response to the removal 
of magnetic energy from the corona. The corresponding frame¬ 
work successfully connects oscillation and contraction, and ex¬ 
plains why some loops oscillate while others do not. 

The key ideas are: 

1. Conversion of magnetic energy reduces magnetic pressure, 
thereby changing the position in which nearby coronal loops 
are in equilibrium. 

2. Since coronal loops possess mass, they cannot respond in¬ 
stantly to changes in their environment. 

3. Comparison of a loop’s natural period of oscillation to the 
time scale of the energy release determines the type of mo¬ 
tion that will result. Loops with periods much longer than the 
time scale of the energy release oscillate around their new 
equilibrium position with an initial amplitude set by the dis¬ 
placement between the old and new equilibrium positions; 
loops with periods much shorter than the energy release time 
scale contract with negligible oscillation; and loops with pe¬ 
riods broadly comparable to the energy release time scale 
exhibit both behaviours at once. 

In conclusion, contraction or displacement of coronal loops 
and transverse loop oscillations are closely related in solar flares, 
and these motions reveal much about the corona and its dynam¬ 
ics, especially rapid decreases in coronal magnetic energy. 
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Appendix A: Coronal loops as harmonic oscillators 

The treatment of coronal loops as harmonic oscillators can be 
rigorously justified for particular equilibria. For example, ne¬ 
glecting gravity, gas pressure, resistivity and viscosity for sim¬ 
plicity and assuming linear perturbations about a static potential 
equilibrium, the momentum equation can be written as 


dfi 


-(V X b) X Bo 

/loPo 


(A.l) 


where ^ is the plasma displacement vector, b is the magnetic 
field perturbation, po is the equilibrium density and Bq is the 
equilibrium magnetic field which we have assumed curl-free. 
Similarly, the time-integrated induction equation is 

b = Vx(fxBo). (A.2) 

We consider a straight equilibrium magnetic field by setting 
Bo = Bqx {Bq is constant to give magnetic pressure balance in 
the equilibrium). Then, performing some algebra, the governing 
equations can be combined to give 


where va = Bq / yjpopQ is the Alfven speed and 'W is a spatial 
differential operator defined by the equation above. In a medium 
with uniform Alfven speed, incompressible solutions (with V ■ 

= 0) to Eq. (A.31 are Alfven waves, while the divergence 
of Eq. ( |A.3[ ) gives the governing equation for fast waves in cold 
plasma. Equation ( |A.3| l is also valid when the Alfven speed is 
a function of position, as happens in equilibria supporting kink 
waves. 

Standing wave solutions are obtained when has separable 
time and spatial dependences, i.e. 


= 7’(f)X(x). 

Substituting this form into Eq. ( |A.3| l, we find 
1 6?-T(t) "WXi(x) 'WX2{x) 


Tit) df2 


Xi(x) 


^2(x) 


= - 0 ) 


(A.4) 


(A.5) 


where X\ and X 2 are the components of X in the chosen co¬ 
ordinate system and we have introduced —aP' as the separation 
constant. It immediately follows that 


dPT 

Af- 


+ oj^T = 0, 


(A.6) 


which is the ordinary differential equation for a harmonic oscil¬ 
lator, while 


WXj(x) ^ -oj^Xjix), je{l,2} 


(A.7) 


implies that Xi and X 2 are eigenfunctions of "W, each having 
eigenvalue -aP. Oscillations appear for real a>, in which case 
Eq. ( |A.6| ) can be scaled so that T represents, e.g., the displace¬ 
ment of the loop apex from its equilibrium position. If u is imag¬ 
inary then the equilibrium can be unstable, but we exclude that 
possibility here for consistency with our observations. The final 
possibility for aP is that it could potentially be complex, which 
may introduce decay alongside oscillation as discussed later. 

The existence of suitable eigenfunctions of "W is expected 
on the basis that standing transverse oscillations of coronal 
loops are observed, including in the event presented in Sect. 

We also provide mathematical justification by considering kink 
modes for the well-known magnetic cylinder model of Edwin &| 
|Roberts| ( 1983| l. In this model, the equilibrium Alfven speed is 
specified in cylindrical coordinates (r, 6, z) as 


ta(i-) 


vao, r < a 

VAe, r > a ’ 


(A.8) 


with VAe > Tao so that the inner region of radius a acts as a 
waveguide for fast waves, giving rise to body modes. Since we 
have Bq constant for equilibrium magnetic-pressure balance, the 


6 


















Russell et al.: Contraction and oscillation of coronal loops 


Alfven speed profile corresponds to a density enhancement in¬ 
side the cylinder. We also assume that there are line-tied bound¬ 
aries at z = 0 and z - L, which gives rise to standing waves with 
wave number k^, the fundamental harmonic having = njL. 
The spatial eigenvalue equations (A.7 1 can be solved for this 


equilibrium using the approach described by Edwin & Roberts 
( |1983| l, which yields Xr{r, 0, z) and Xg(r, 0, z) with radial depen¬ 
dences in terms of Bessel functions. Meanwhile, the eigenvalue 
is determined by 

(A.9) 


O) 


- r^k^ 


where Cp(ak^) is the phase-speed of the kink wave. In the thin- 
tube limit, k^a <si 1, the waves become non-dispersive and Cp 
reduces to the kink speed. 


Ck 


Povlo +PevlA 

PO+Pe j 


1/2 


(A.IO) 


(A.ll) 


d^x 


dx 


—zir -I- An {x - 5^0(0) + = 0, 


dP- 


dt 


(B.l) 


where t — tjT with r = Inju), and x - xjD with D an appro¬ 
priate length scale such as the final displacement. This equation 
is equivalent to a pair of coupled first-order ordinary differential 
equations. 


dx 

dt 

dv , 

— = -Airix - xo(t)) - AnKV, 
dt 


(B.2) 


which are easily integrated using fourth-order Runge-Kutta. 
Setting the damping parameter as /c = 0.1 produces decay simi¬ 
lar to the motivating example. Solutions of the three types shown 
in Fig.fflcan then be obtained using different driving functions 
for xo(fk Impulsively excited oscillations (Fig. [^) result if the 
system is driven using a step function, e.g. 


Since spatial eigenfunctions have been found and up is real and 
positive, standing kink waves in the classic magnetic cylinder 
model are consistent with a harmonic oscillator view of coronal 
loops. 

More generally, standing kink oscillations can decay over 
time due to physical damping by viscosity, resistivity etc. or 
due to transfer of energy from the kink mode by resonant ab¬ 
sorption ( Ruderman & Roberts |2002| Goossens et al.|2002) l. In 
these situations, the above approach of separation of variables 
can still be applied, leading to an oscillation equation of the form 
of Eq. ( |A.6[ ). However, solution of the spatial operator equations 
now yields a complex frequency, co - u)r + ia>i. When this occurs 
it can be useful to recast the oscillation equation as a damped 
harmonic oscillator with real coefficients; 


xo(f) = 


2, f < 4 

1, t > tc 


(B.3) 


where P is the time of the collapse. Examples of gradual dis¬ 
placement (Fig.|^) can be produced using 




3 - tanh 


t-tr 


(BA) 


with A > 1. In fact, Eq. (B.4i can produce any of the three 
types of response, depending on the value of A, and converges to 
Eq.Q for A ^ 0. Nonetheless, the clearest examples of os¬ 
cillation during collapse (Fig.|^) are produced when the change 
in equilibrium starts sharply, e.g. under 


dPT , dr 

— +co,T+2ujok— - 0, 

where ojq is the frequency of corresponding undamped oscillator 
and K is the damping ratio. To see the correspondence, observe 
that harmonic solutions of this new equation with an e'"' time 
dependence also have complex frequency, with atr - wq Vl - 
and a>j - ujqk (we will consider underdamped solutions with 
0 < a: < 1 so Wr is always real and cu, > 0 giving decay). This 
implies that the original oscillation equation with complex at has 
the same oscillatory and decay behaviours as the damped oscil¬ 
lation equation with real coefficients set as Wq = wj -H and 
a>oK - u)i, where cOr and w, are obtained from the spatial eigen¬ 
value equations. 

We conclude by noting that the approach described above 
is fairly general and it can be applied to analyses that make 
fewer simplifying assumptions. For example, gas pressure was 
in fact included by Edwin & Roberts ( |1983| l, while magnetic 
curvature, density stratification, flux tube expansion and non¬ 
circular cross-sections have all been considered since then (see 
Ruderman||20031 |V an Doo rssela ere et al.| 20041 [Andries et al.| 
2005 1 Ruderman et al.|2008||Ruderman|2015| l. In each case, suit¬ 
able spatial eigenfunctions were obtained and a dispersion rela¬ 
tion giving u) was derived, which determines the real coefficients 
in the corresponding damped oscillation equation. 


Appendix B: Solution of the model equation 

Our model equation produces all three types of solution dis¬ 
cussed in Sect.[^ Writing Eq. Q in dimensionless form. 


xo(t) = 


2, t < tc 

2 - tanh ((f - fc) /A), t > tc 


(B.5) 


with A = 1. 


Appendix C: Effect of switch-on time 

It was noted in Sect. that the greatest amplitude oscillations 
are obtained when the change in equilibrium position is ini¬ 
tially sharp. Physical justification was given in Sect. but a 
quantitative demonstration is also provided here. Using the non- 
dimensionalised model of Eq. (B.l I with damping turned off 
(k - 0), we consider a test problem in which the equilibrium 
position accelerates from rest to motion at a constant speed. 
Assuming the equilibrium position accelerates at a constant rate 
during a switch-on interval equivalent to 6 periods of oscillation. 


xo(f) = 


Xr-fl( 26 ), tc <t <tc +5 . 
Xr — t + (5/2, tc + 6 <t 


(C.i) 


Here, the normalising length has been set to the distance that 
the equilibrium position moves per period in the constant speed 
phase, t > tc + 6 . The specified driver produces oscillations su¬ 
perimposed on a never-ending collapse, and the amplitude of 
the oscillations is readily measured from the long term solution. 
Figure [CT| plots the normalised amplitude for a range of switch- 
on times, 6, which is the only free parameter in the dimensionless 
model. In this test problem, dimensional amplitude is limited 
only by the maximum rate of collapse, with amplitudes capped 
at 0.16 times the maximum displacement per period. Nodes in 
Fig. C.I show cases where acceleration of the equilibrium posi¬ 
tion resonates with the loop period in such a way that long-term 
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Fig. C.l. Amplitude of oscillation as a function of the energy- 
release switch-on time, for a test problem in which the equilib¬ 
rium position accelerates from rest to motion at a constant speed. 
The switch-on time (duration of the acceleration) is normalised 
to the period of oscillation, and the amplitude of oscillation is 
normalised to the equilibrium displacement per period after the 
acceleration. 


oscillations are not produced. Most significantly, the plot con¬ 
firms that the largest amplitude oscillations are excited when the 
switch-on time is a small fraction of the oscillation period. 

Examining the motion of loop L3 in Fig. the inferred 
equilibrium position shifts by approximately 12 Mm during the 
loop’s first period of oscillation. Multiplying this distance by 
0.16 gives an estimate of 2 Mm for the maximum amplitude that 
can be excited by the corresponding rate of contraction. That es¬ 
timate is very close to the actual amplitude seen for L3. Based on 
these values and the fall-off of amplitude with increasing switch- 
on time (Fig. C.l i we suggest that a quarter of L3’s period (ap¬ 
proximately 40 s) is a reasonable upper limit to place on the 
energy-release switch-on time for the SOL2012-03-09 flare. 
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